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Abstract 

The homogenization of a metamaterial made of a collection of scatterers pe- 
riodically disposed is studied from three different points of view. Specifically 
tools for multiple scattering theory, functional analysis, differential geome- 
try and optimization are used. Detailed numerical results are given and the 
connections between the different approaches are enlightened. 
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1. Introduction 

Giving a general definition of what homogenization is is difficult, because 
of the various meanings attached to it. The physics at stake is the situa- 
tion when a wave illuminates a complicated object, generally consisting of 
a periodic set of scatterers, contained in some domain Q and gives rise to a 
diffracted field U s . Loosely speaking, the homogenization problem consists 
in identifying homogeneous constitutive relations such that the same domain 
fl inside which these constitutive relations hold, leads to a diffracted field 
such that U s and are close to each other, in some meaning to be speci- 
fied. This can be done reasonably only if the wavelength A is larger than the 
period d. This identifies a small parameter r\ — d/X. 

Having been educated by mathematicians, the definition that I would 
consider the best one is the following: consider a partial differential equation 
P n with oscillating coefficients a n {x). Consider a solution u v of the equation: 
P v {u v ) = f, where / is some convenient source term. Then the goal of 
homogenization theory is to find a convenient topology in which u v converges 
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to a function u satisfying an equation: Po(u ) = f. The operator P is called 
the homogenized operator This definition is quite clear and at the end, 
it leads to results such as: "when 77 tends to 0, tends to uo in some 
specific meaning" . The point of being able to specify a convergence is very 
interesting, in that it gives a clear meaning to the question "how close are 
f/ s and f/* ?". 

In the metamaterials community, it is not rare that the notion of parame- 
ters extraction be used as a homogenization scheme [2J. In that situation the 
structure is considered a black box (or rather a black slab!). This pragmatic 
approach, although it might be useful, cannot be considered a homogeniza- 
tion procedure. 

Sometimes the Bloch spectrum is used as well. However, care should be 
taken because of the following result: 

Proposition 1. Given any isofrequency dispersion curves given implicitly 
in the form F(k x , k y ) = 0, there exists a spatially and temporally dispersive 
permittivity e(k,u>) reproducing these curves. 

Proof. Write: k 2 = k 2 + F(k x ,k y ) and define: e(k, uj) = k +F j^ x ' k y^ ; where 
k = u/c. □ 

This shows that given any dispersion curves, it is always possible to re- 
produce it by using a spatially dispersive permittivity. However, nothing can 
be said on the complete electromagnetic field inside the medium: the Bloch 
diagram only accounts for the plane wave part of the field. The introduction 
of spatial dispersion should only be made with great care 0] . 

In the following, the structure considered as a model problem is a periodic 
set of 2D scatterers (cf. fig. (pQ)). The medium is infinite in the x direction 
and it is made out of a stack of basic layers made of an infinite number of 
rods periodically disposed at points x = nd. 

In the following, three different points of view are given on the homoge- 
nization of this structure: multiple scattering, double-scale, micro-local. 

2. Multiple scattering homogenization 

The approach proposed here is reminiscent of that in [4|, although it is 
rather different because here no averaged field is defined. It would be inter- 
esting to make a connection between these two approaches. 
Our point is to describe the electromagnetic behavior of one line of scatterer 
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Figure 1: Sketch of the structure under study. It is made of a stack of gratings, each 
consisting of periodically disposed scatterers. 



(a grating) alone. Each scatterer "n" is characterized by a scattering coeffi- 
cient s . When it is illuminated by an incident field E\ it gives rise to a field 
u s (r) = bH (ko\r\) where b = s E l (n). For the infinite set of scatterers, this 
gives a diffracted field that reads: 



u s 



J2KH (k \r-r n \) (1) 



the field diffracted by the n th scatterer is obtained by writing that it is the 
response to the incident field E l and the field difffracted by the other scat- 
terers: 



b n = s I E\n) + ^ b m H (k \m - n\d) ) 



(2) 



If the incident field is pseudo-periodic, i.e. E l (n) = e lkd E % (n — 1), then it 
holds: 

b = (3) 

° ^-J2 m ^e^H (k \m\d) 
the series that enters this relation can be written: 

£ e^H (k \m\d) = -1 - + * In (£) + A- 

it it \k dj dp 
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An asymptotic analysis of this series [5[ allows to write the following expan- 
sion : 

Lemma 1. 

^(rt-4H (M)2) (4) 

Above the grating, the propagative part of the electric field reads as: 
E+{y) = e~ iPoy + r(>o, A))e i/3oy and below it reads E~(y) = t(k , l3 )e- il3 ° y 
where: r(k Q , (3 ) = and t = l + r 

Energy conservation implies that : |r| 2 + |1 + r| 2 = 1 and therefore the 
following representation holds: 

Proposition 2. There exists a real function xi^o, Po) such that: 



l + ix(k ,(3 ) 

Proof, energy conservation shows that: = —1, the theorem follows by 

defining: x = w- D 

Proposition 3. For the propagative part of the field, the grating is equivalent 
to an infinitely thin slab whose transfer matrix is: 

where L = x/2/3 - When the ratio a/A is very small, it holds: L ~ ^ In (^-) 

Proof. The electric field is continuous, the matrix T g is obtained by comput- 
ing the jump of the normal derivative of the field: d y E + (0)— d y E~(0) = lifter. 
Using proposition 1, we get: E(0) = 1 + r = = —ix r - Hence: 

2if3 r = -^E(0). ' □ 

A layer of the structure can thus be described by this slab surrounded by 
homogeneous layers of height h/2, the transfer matrix of this basic sandwich 
structure (cf. fig. ©) being: T = T(h/2)T g T(h/2), where : 

T (h/9\-f C0S (/W 2 ) sm(/3 h/2)//3 \ , , 

1WZ) I -A>sin(A,/i/2) oos(A)V2) J [ ) 
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Figure 2: The homogenization of the sandwich structure. The horizontal black line 
represents the grating. 



A simple calculation shows that: 



O sin(/3 h) co S ((3 h) + s -^ 

The transfer matrix is known to be an unstable object for numerical pur- 
poses. The preferred quantity is the scattering matrix j6[. Therefore, the 
homogenization is now performed numerically by searching for a permittiv- 
ity e e ff(k , (3 ) such that the scattering matrix of a homogeneous slab with 
this permittivity fits that of the sandwich structure. Specifically, the cost 
function is: 

J[e c s] = \r g - r eS \ 2 + \t g - t cS \ 2 (8) 

where r g and t g are the reflection and transmission coefficients of the sandwich 
structure and r eS ,t eS that of the homogeneous slab. 

2.1. Numerical applications 

First the scatterer in the basic cell is a dielectric rod with permittivity 
e r = 9 and radius a/d = 1/2, i.e. the rods are touching. In that case, the 
low-frequency behavior is well-known: the sandwich structure is equivalent 
to a homogeneous slab with permittivity: £i w = 1 + (e r — l)ir (^) . 
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Figure 3: The reflection spectrum for the homogenized structure (dashed) and the sand- 
wich structure (continuous) . 

This serves to test the proposed approach and also to evaluate the spatial 
dispersion effect in the medium. The reflection spectra for a plane wave in 
normal incidence and for both the sandwich structure and the homogeneous 
slab are given in fig. ([3]). It can be shown that the homogenization works 
very well for X/d > 15. The resulting homogenized permittivity, depending 
on both the frequency and the horizontal Bloch vector, is given in fig. (Tj0). The 
averaged value is £i w = 7.28318 ± 1CT 5 . The value obtained numerically for 
X/d = 1000 is: e e ff = 7.2832. The homogenized permittivity was calculated 
for two angles of incidence and 7r/4: it can be seen that the effects of spatial 
dispersion are quite small. 

Second, we choose a resonant scatterer with a scattering coefficient of the 
form 

1 u - u z 

so = , (9) 

where g(u) is regular and satisfy g — > as to — > 0. This pole-and-zero 
form is quite a common one. The sandwich structure air-grating-air can be 
replaced by a slab with permittivity e e g(u, (3o) obtained from the optimisation 
procedure described above. The reflection spectrum is given in fig. ([5]), the 
homogenization procedure works very well for X/d > 12. The corresponding 
homogenized permittivity is given in fig.(j6]). The spatial dispersion effects are 
negligible in that situation. Interestingly, it seems that resonant structures 
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Figure 4: The effective permittivity obtained by using the cost function ([8]). The con- 
tinuous curve corresponds to /3o = and the dashed one to /3o = ko^/2/2. The horizontal 
line corresponds to the averaged permittivity £i w . 




Figure 5: The effective permittivity obtained by using the cost function (j8]) for a grating 
of resonant scatterers with scattering coefficient given by ([9]). 
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Figure 6: The effective permittivity obtained by using the cost function (JH| . 



can be homogenized at smaller ratio X/d than non resonant ones. 



3. Pseudo-differential homogenization 

In this section, a homogenization approach inspired by j3] is given. For 
the sake of simplicity, the theory is specialized to the Helmholtz equation: 
Au + k^eu = f, where e is the dielectric function of the structure in fig.([TJ). 
Let us define by T the operator such that: u = T(f). It is a pseudo- 
differential operator [8] whose symbol is denoted a(r, k). The following form 
of u holds: 

u(r) = [ a(r,k)f(k)e ih - r dk (10) 



The symbol is the response of the system for a plane wave source. It can be 
shown that it has an expansion in terms of the Bloch spectrum 

a(r k) y \Mr,k)\ 2 

a ^~l.E p {k)-kl { ] 



where (f) p is such that: <p p e r is periodic and — s~ A(f) p = E, 
Let us now consider the periodic medium with basic cell Y . The polariza- 
tion field is: = Eo(e — l)u v . Define the average field through: u e s = 
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J Y a 7?( r ; k)e lk ' r dr and the average of the polarization field: 

P c s = £o / (e — l)s( r > k)e tkr dr = e (e - 1) / ^(r, k)e lk ' r dr 
Jy Jd 

The effective permittivity is then £ eff (fc, u;) = 1 + £ ^ ff , that is: 

^(*,^)=i+( e -i) {°^ r '*|Cf < 12 ) 

Jy o, v (r, k)e ZK r dr 

It is interesting to note that the expression above is regular whatever k, 
despite the fact that the symbol has poles. In the very low frequency domain, 
one obviously has: 

eeff(A;,a;) = l + (e- (13) 

This approach could probably be made better by using a multiple scale 
expansion of the symbol in the form: 

a n {r, k) = a (r, r/rj, k) + i qa 1 {r,r/r], k) + ... 

. Besides, it would also be interesting to recover the classical homogenization 
results for the case of a magnetic field linearly polarized along the wires. Work 
is in progress in that direction. 

4. Multiple scale homogenization 



This final approach is largely used in the mathematical community [15 . 



16] . also it is practically ignored by physicists [10J. By its very definition, it 
is the only one that can give a clear meaning to the notion of convergence of 
the fields. I will only give a sketch in a well-known situation, but I will give 
a new framework by using differential forms [l[ , because the structure of the 
theory is very nice then. Maxwell equations are written by using differential 
forms: 

dE v = iojfMoH* 1 , dH v = -iue v * E v (14) 

where d denotes the exterior derivative, E v = E^dx n and H n = H%dx n , n = 
1,2,3 are 1-forms (obtained from the usual vector fields by means of the b 
operator and a flat metric and * stands for the Hodge operator. The 
starting point is to assume the following expansions for E v and H v : 

E v (x) = E°(x, x/rf) + t] E 1 (x, x/rj) + r] 2 E 2 (x, x/rj) + ... , , 
H v (x) = H°(x, x/rj) + r] H x {x, x/rj) + r] 2 H 2 (x, x/rj) + ... 
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where the fields E°, H°, E 1 , H 1 ... depend upon two sets of variables (x,y) e 
R 3 x Y. 

This can be justified in the framework of double-scale convergence (lil . lil 
but this would take us to far. 

There is a nice mathematical structure linked to the multiscale expansion. 
The projection: 

ir:E (x,y) — y J E {x,y)dy (16) 

defines a fiber bundle. Where the fiber is the cotangent bundle to the flat 
torus: F = T*(M 2 /Z 2 ) and the base is the cotangent bundle to the ambient 
space :B = T*(M 2 ).Its trivialization is as follows: 

7r -1 (Z7) ~ FxB (17) 

where U is an open neighborhood in the base. The y variable plays the role 
of a hidden variable that accounts for the microscopic behavior inside the 
basic cell. 

The expansion ( IT5|) induces the following splitting of operator d: d = 
d x + ^d y . We arrive directly at the systems satisfied by the microscopic 
fields, holding on Y: 

d y E = J d y B = 



5 y (eE ) = \ 5 y B = (18) 

where 5 = — * d* is the exterior co-derivative. The first system involves 
the electric field alone and is of purely electrostatic nature. Because the 
limit field only depends microscopically on (2/1,2/2), one immediately obtains: 
E®(x,y) = Es(x). Moreover the transverse microscopic electric field E\ = 
{Ei, E%) is exact on the basic cell Y. This implies that Ej_{x, ■) reads as: 

E° ± {x,-) = E ± {x) + d y ip{-) (19) 

and therefore the effective electric field E±{x) = {Ei{x) , E2{x)) belongs to 
the first de Rham co homology space H 1 {Y), a space isomorphic to IR 2 . The 
complete solving of the miscropic system is done by a linear decomposition: 
ip = Ei<fi + E 2 (p2, where the functions cpj,j = 1,2 satisfy: 

y, = 0onY\D 
Un^j — ~ n j on oD 
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This leads to the linear relation: E = S(y)E(x), where: 




£(y) = d 2 <pi 1 + d 2 <p 2 (21) 



The second system shows that the magnetic field is both exact and co- 
exact, implying that B does not depend on y. We denote: B (x,y) = B(x). 
Finally, the macroscopic equations read as: 

dE° = iunoH , dH° = -iue eE° (22) 

After averaging on Y, one obtains: 

dE = iufi H, dH = —iueoe e sE (23) 

where e c q = f Y e(y)S(y) dy. The usual anisotropic permittivity tensor is 



found. Generalizations and details can be found in [171. |18|. The method can 
be extended to deal with resonant structures Jl9j] and obtain a homogeniza- 
tion result for higher bands that the first one, also it is sometimes wrongly 
believed that this approach can only deal with quasistatic problems. It is 
interesting to note that the degree of the forms involved can be a clue to the 
definition of the averaged field 2^ . 



5. Conclusion 

We have described three different approaches to the homogenization of 
a two dimensional dielectric metamaterial. The micro-local approach is less 
developed than the other but still seems quite interesting, and in need of 
mathematical development. Our preferred one still is the multiple scale ap- 
proach because it can deal with the notion of convergence and the boundary 
conditions naturally. However, it is also in need of mathematical development 
in order to take spatial dispersion into account. 
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